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INVESTIGATION OF AN ANOMALOUS FLOW CONDITION OF THE 
LANGLEY PILOT MODEL EXPANSION TUBE 


By Wilfred J. Friesen 
Langley Research Center 

SUMMARY 

Free -stream flow velocity measurements were made in the Langley pilot model 
expansion tube during the test flow interval. During this interval, an anomalous dip in 
pitot pressure occurs for the expansion tube operating conditions employed. Within the 
test flow interval, the main conclusions reached from comparison of the measured flow 
velocity, pitot pressure, and tube wall pressure are: the variations which occur in veloc- 
ity and wall pressure are small compared with the variations in pitot pressure; a corre- 
sponding dip in the derived flow density is associated with the dip in pitot pressure; and 
the value of the average density over the interval, which results from the expansion from 
the shocked intermediate chamber condition, is approximately one-half of the value that 
can result from only an isentropic process. 

INTRODUCTION 

An expansion tube is an impulse device which produces a short- duration, hypersonic 
test flow with flow properties suitable for the study of problems related to planetary entry. 
The results of an exploratory study of the performance of the Langley pilot model expan- 
sion tube are reported in reference 1. 

In reference 1, an unexpected decrease or dip in the pitot pressure was reported. 
The pitot pressure dip occurs within the test flow time interval for certain operating 
conditions of the expansion tube and seriously reduces the normally short testing time of 
the facility. Interpretations of this dip in pitot pressure depended on the assumptions 
made concerning the test flow velocity and density. The velocity of the interface (which 
is the first flow boundary of the test flow to arrive at the test location) was measured, but 
measurements of flow velocity within the subsequent test flow were not available to Jones 
and Moore. (See ref. 1.) 

The purpose of the present investigation was to derive the flow density associated 
with an anomalous flow condition from measured values of the test flow velocity on the 
tube center line, pitot pressure, and tube wall pressure. Although specifically directed 
toward a better understanding of the nature of the anomalous flow, this investigation is 



part of an effort directed toward understanding the cause of the anomaly. The scope of 
the present investigation was limited to one initial operating condition of the expansion 
tube for which a significant pitot pressure dip occurred within the nominal test flow 
time interval; and only nitrogen was used as the test medium (similar pitot pressure dips 
are observed for other test gases as well). Also a necessary part of the present investi- 
gation was the development of a method of measuring the flow velocity. 

The results of the measurements of flow velocity, pitot pressure, tube wall pres- 
sure, interface velocity, and shock velocity that were obtained for a specific anomalous 
flow condition of the expansion tube are presented. Also presented are results which 
were derived from the measured values. These results include the density, temperature, 
integrated flux, and the reproducibility of several of the flow properties of the expansion 
tube cycle. 

The magnitudes of the flow properties of an expansion tube can be varied over large 
ranges. To aid the reader who may be unfamiliar with the nominal magnitudes and trends 
encountered in expansion tube flows, the predictions of a simple, highly idealized flow 
model are presented. 


SYMBOLS 


a speed of sound 

I intensity 

Iq initial intensity 

M flow Mach number 

n number of expansion tube runs 

p static pressure 

pj. pitot pressure 

P w tube wall pressure 

T temperature 

t time 
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* 0 , 11 * 0,2 

times corresponding to positions 1 and 2 

^*c = *c,2 

" tc,l 

*P 

time corresponding to peak intensity 

Us,l 

velocity of primary (nitrogen) shock wave 

u 

flow velocity 

u i 

interface velocity 

u i,P 

average interface velocity for pressure runs 

um 

measured flow velocity 

X 

h-L 

to 

positions 1 and 2 (slit centers) 

Ax = x 2 - 

X 1 

y 

specific heat ratio 

p 

density 

a 

standard deviation from the mean 

a 

estimated variance of the mean 


Subscripts: 

1,2,5, 10, 20 refer to flow regions indicated in figure 1 

IDEALIZED MODEL OF FLOW CYCLE 

A definition and description of an expansion tube, its operating flow cycle, theoret- 
ical basis and performance capability are discussed in reference 2. 

In this section an idealized model of an expansion tube flow cycle is briefly sketched 
and some of the predicted flow properties are stated. This one -dimensional model 
assumes perfect gases, nonviscous flow, isentropic expansions, no mixing of gases, and 
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idealized diaphragms. The predictions of this model are presented primarily as a guide 
to indicate only the expected magnitudes and trends for an ideally operating expansion 
tube. 


Flow Cycle 

A distance-time diagram of an idealized flow model is shown in figure 1. Initially, 
three gases are at rest in three chambers, which are defined by two diaphragms in a 
constant area tube. The flow cycle is initiated by the rupture of the first diaphragm. 

After the rupture of the first diaphragm, the flow cycle, up to the second diaphragm, 
is that of an ideal shock tube. The expansion of the hydrogen driver gas results in the 
compression and acceleration of nitrogen within the intermediate chamber at a constant 
supersonic rate. The flow velocity and state properties of the nitrogen throughout 
region (2) are constant. Region (2) is bounded by the nitrogen shock and the nitrogen- 
hydrogen interface (contact surface). The secondary diaphragm ruptures upon arrival of 
the shock wave. 

After the rupture of the second diaphragm, the velocity of the nitrogen-hydrogen 
interface remains constant as the moving nitrogen expands into the acceleration chamber. 
For a brief time the flow cycle in the acceleration chamber can be viewed as a shock 
tube cycle with the nitrogen "driver" initially in motion. Since the motion of the nitrogen 
is supersonic, the upstream facing expansion wave which moves into region (2) is washed 
downstream. The test gas is that portion of the nitrogen which has expanded into 
region (|). Throughout this region the nitrogen flow velocity and state properties are 
constant. Region (§) is bounded by the helium -nitrogen interface (contact surface) and the 
tail of the expansion fan. 

The role of the helium acceleration gas is to provide a boundary which moves at 
constant velocity and limits the expansion of the nitrogen to a finite constant pressure 
P 5 < The conditions of constant velocity and pressure require that the work done on the 
boundary should proceed at a constant rate. For this idealized model, the necessary 
boundary work is accomplished by continual compression and acceleration of helium at a 
supersonic rate using the nitrogen test gas (region (5)) and the compressed helium 
(region as the contacting or pressure transmitting media. 

Flow Properties 

The following sequence of flow properties is predicted for an observer located at 
the test section. Flow begins with the arrival of the helium shock. A helium flow with 
constant velocity and constant state properties is observed until all the helium, which was 
initially located between the second diaphragm and the test section, has passed the test 
section; the nitrogen test flow with constant velocity and constant state properties is now 
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observed until the test flow is terminated by the arrival of the tail of the expansion fan. 
After the arrival of the expansion fan, the flow properties of the nitrogen are no longer 
constant with the passage of time. During the passage of the expansion fan the flow veloc- 
ity decreases and the state properties increase monotonically. 

In general, the following relations between the helium flow and the nitrogen test 
flow are observed at the test section: 


u 20 = u i = u 5 = Constant 
P 20 = Pg = Constant 

p 20 * p 5 
t 20 * T 5 


Since the expansion of the nitrogen from region (2) to region (fT) is assumed to be an 
isentropic process, 


£5 = (*z\ 1/Y 

P2 \ p 2 / 




( 1 ) 


The relation between the flow velocities in region (2) and region (T) for the one 
upstream facing expansion wave, which is given in reference 2, is 


y 


Tf a 2 + u 2 



a 5 +u 5 


( 2 ) 


APPARATUS 
Expansion Tube 

The Langley pilot model expansion tube has been completely described in refer- 
ence 1. Only one set of chamber lengths and one set of gases and initial pressures were 
used in this investigation. The gases used in the driver, in the intermediate chamber, and 
in the acceleration chamber were, respectively, room temperature hydrogen, nitrogen, 
and helium. The first diaphragm was a scribed, 1.59-mm-thick steel diaphragm and 
the second, a 0,0064-mm-thick mylar diaphragm. At the test section the flow exited as 
a free jet and subsequently entered a large dump tank. The dimensions and initial 
operating conditions for the tube are given in table I. 
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Measurement Apparatus 


Tube instrumentation .- The expansion tube instrumentation that was used to detect 
wave passage and provide trigger signals has been described in reference 1. Ion gaps 
and pressure gages were used to detect the nitrogen shock in the intermediate chamber; 
pressure gages were used to detect the helium shock in the acceleration chamber. A 
microwave antenna located downstream of the test section was used to track the helium- 
nitrogen interface in the acceleration chamber. A photodetector was used to detect the 
arrival of the helium -nitrogen interface at the test section. 

Improvements to the facility vacuum system, instrumentation, and operating pro- 
cedure have occurred since the exploratory work reported in reference 1, and have 
resulted in improved measurement of the initial helium pressure and a substantial reduc- 
tion in the contamination level of the helium. With regard to contamination, for example, 
no light was observed from stagnated regions of the helium flow during any of the tube 
runs. Use of an improved absolute pressure gage resulted in a reproducible filling of the 
acceleration chamber with the low pressure helium; the pressure gage was of the capac- 
itance diaphragm type. 

Flow pressure .- Piezoelectric pressure transducers were used in the measurements 
of tube wall pressure and pitot pressure. The wall pressure gage was located approxi- 
mately 26 cm upstream of the test section and was mounted flush with the tube wall. The 
unfiltered signal from this transducer was recorded by an oscilloscope. The pitot pres- 
sure probe was located at the test section on the tube center line. The probe, which had 
an external diameter of 0.95 cm, shielded the transducer from particles in the flow. The 
signal from this transducer was processed through a passive filter to suppress frequen- 
cies above 50 kHz and then recorded by an oscilloscope. The locations of the pressure 
probes relative to the exit of the expansion tube are shown in figure 2. 

Flow velocity .- A necessary part of this investigation was the development of a 
method of measuring the flow velocity suitable for the flow conditions to be encountered. 
The method used consisted of producing a luminous region within the flowing gas, and 
then, by means of photodetectors, observing the time of arrival of the region at two 
fixed downstream positions in the flow field. This method is similar to those reported in 
references 3 and 4. The specific differences are primarily involved with the techniques 
for producing the luminous regions. In the present investigation a localized region in the 
flow was excited to luminescence by passing a short duration (10" ^ sec) current pulse 
through a previously ionized column in the gas. The ionized column was formed by 
photo ionization just prior to its excitation. Photoionization was used to identify or define 
a region within the flowing gas for the velocity measurements reported in reference 5; 
however, the method employed to detect the subsequent location of the defined region was 
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not applicable to the larger flow densities encountered in the present investigation. Since, 
for these larger densities, the ionization produced decreases rapidly with distance from 
the photon source, it was necessary to employ a spacing of about 2.5 cm between the 
photon source and the exciting electrode and to design both of these components to operate 
within the boundary of the flow field. 

A schematic drawing of the windowless ultraviolet light source used to photoionize 
the column is shown in figure 3. Ultraviolet light (UV) was produced by means of a con- 
stricted electrical discharge in the nylon capillary which was 0.78 mm in diameter and 
1.6 cm long. The light duration was of the order of 10 ^ second. Light from the capil- 
lary was partially collimated by a 0.78-mm-diameter field stop. The collimator assembly 
also served as the electrostatic shielding and was operated at ground potential. The 
exciting electrode consisted of a 0.5-mm-diameter tungsten wire. The distance between 
the exciting electrode and the exit collimator of the light source was 2.2 cm. The light 
source and electrode were located in a plane through the expansion tube center line, and 
their locations relative to the exit of the expansion tube are shown in figure 2. 

A circuit diagram of the transmission line pulse generator which provided pulses 
to both the light source and exciting electrode is shown in figure 4. The generator pro- 
duced a pulse of about 10 ^ second duration. About 4 joules were stored in the coaxial 
transmission lines when charged to 40 000 volts. The two pulses were derived from the 
same pulse generator in order to insure proper synchronization between photoionization 
and excitation of the gas column. The pulse to the exciting electrode was delayed about 
10" r7 second with respect to the pulse to the light source by means of an additional length 
of coaxial cable. 

The passage of the luminous column at two fixed positions in the flow field was 
detected by two photodetectors (type 931- A); these photodetectors were located behind 
rectangular slits. The optical arrangement used is shown in figure 5. This optical 
arrangement effectively placed two rectangular slits (1.59 mm wide, 3.13 mm long, and 
1.905 cm apart) in a plane through the center line of the expansion tube. The upstream 
slit was located 1.3 cm downstream of the photoionizing source. The lens was used at a 
magnification of two and had a focal length of 17.8 cm and an aperture stop of f/2.5. 

An optical filter was used to block radiation to the photodetectors in the wavelength 
region below 5000 A, in order to remove radiation due to the first negative band system 
of the excited N* ion. Radiation from this band system normally accounts for a large 
fraction of the light emitted from an electrical discharge in nitrogen. From the stand- 
points of photodetector overload requirements and interpretation of the data, the com- 
bined high intensity and rapid decay associated with this radiation was undesirable for 
the present application. 
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The photodetector outputs were terminated in 50-ohm load resistors, and their sig- 
nals were recorded on two channels of an oscilloscope. An overall rise time of about 

_7 

10 second is estimated for the photodetector system. Within the range of light intensi- 
ties encountered, the responses of the photodetectors were nonlinear to changes in light 
intensity. 


DATA ANALYSIS 

The measured results presented in this paper represent average values obtained 
from 30 expansion tube runs for a single operating condition. Of these, 4 runs were used 
to obtain wall pressure and pitot pressure data, and 26 runs were used to obtain flow 
velocity data. A tube run was required for each determination of flow velocity. The flow 
velocity measurements were made during the nitrogen flow at seven nominal times fol- 
lowing the passage of the helium-nitrogen interface. At least two expansion tube runs 
were made for each nominal time; the reproducibility of the nominal time locations was 
about ±5 psec. 

Since small variations occurred in the opening of the primary diaphragm, a means 
of comparing the flow velocities obtained from different expansion tube rims was required. 
The variation in the measured helium-nitrogen interface velocity Uj was used to char- 
acterize the resultant variation in the measured flow velocity u m due to variation in 
the primary diaphragm opening. Each measured velocity determination u m was nor- 
malized by the corresponding measured interface velocity ui, and then multiplied by the 
average interface velocity u^p obtained for the four pressure runs. The average 
velocity u which is presented for each of the time locations was obtained from 


The standard deviation from mean interface velocity for the four pressure runs was 
±0.8 percent. The standard deviation for the 30 runs was ±2 percent. No normalization 
correction was applied to the pressure data. 

Pressure 

Sample data oscillograms of the tube wall pressure and pitot pressure are shown in 
figure 6. The data were reduced to numerical form for computation at time increments 
of 5 p sec. 

Calibration of gages .- The two pressure gages were calibrated by relating the initial 
rise of the transducer signals to the corresponding pressure changes associated with the 
passage of a normal helium shock in the acceleration chamber. The pressure changes 
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were calculated by use of the adiabatic relations for a normal shock given in reference 6. 
The required inputs to the calculations were the initial conditions of the helium and the 
measured helium shock velocity. Helium was assumed to behave as a perfect gas. 
Because of the finite rise time of the transducer signals, the initial signal values fol- 
lowing shock arrival were obtained by linearly extrapolating the trend of the signal data 
to the time of shock arrival (t = 0). For the extrapolations, data were used in the range 
of 5 to 55 iLisec for the wall pressure gage and 50 to 150 /isec -for the pitot pressure gage. 
The standard deviations from the mean values obtained for the four expansion tube runs 
were; ±1 percent for the two pressure changes, ±1.5 percent for the wall pressure cali- 
bration factor, and ±3 percent for the pitot pressure calibration factor. 

Possible systematic errors.- The largest uncertainties in the reported pitot and 
wall pressure measurements are expected to be of a systematic nature. Uncertainties in 
the degree of validity of at least two assumptions required in the gage calibration, lead to 
possible systematic errors of the order of 10 percent. It is possible that the values 
reported for the tube wall pressure and pitot pressure associated with the nitrogen flow 
may be as much as 10 percent larger than the true pressures. 

It was assumed that the pressure response of the gages was linear over the entire 
range of pressure levels encountered by the gages. A nonlinear behavior would be 
expected to introduce the larger error in the pitot pressure measurements, since this 
probe encounters a much larger change in pressure level because of the change from 
helium to nitrogen flow. The pressure levels associated with the helium flow were below 
the pressure range for which the linearity of gages of this type had been tested. A pos- 
sible systematic error in a range from 0 to 10 percent for pitot pressure measurements 
during the nitrogen flow was estimated from extrapolation of linearity data for gages of 
the same type. 

It was assumed that the tube wall pressure was equal to the free-stream static 
pressure on the tube center line. The degree to which this assumption is valid is not 
known. The degree of validity may also change for flow immediately following the helium 
shock and subsequent helium flow due to boundary-layer growth. 

Based on a static calibration of the transducer, the initial rise in the observed tube 
wall pressure was about 12 percent smaller than the free-stream pressure calculated 
from the helium shock velocity. The discrepancy could be reduced to about 8 percent by 
considering a probable attenuation of the helium shock in the acceleration chamber of the 
order of 2 percent. It was not possible to determine whether the remaining discrepancy 
resulted from a nonlinear behavior of the transducer or a real aerodynamic effect. (The 
pitot pressure transducer was statically calibrated but an error in calibration was dis- 
covered too late to recalibrate this transducer.) 
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Based on reference 7, an error of less than 3 percent is expected to be introduced 

in the pitot pressure measurements because of the change in Reynolds number as the flow 

changes from helium to nitrogen. Estimates of the unit Reynolds numbers are 5 x loVm 

fi / 

for the helium flow and 2 x lOym for the nitrogen flow. 

Shock and Interface Velocity 

Intermediate chamber shock velocity .- The reported nitrogen shock velocity is the 
average velocity over the distance interval located approximately from 2.1 m to 0.4 m 
upstream of the secondary diaphragm. The estimated precision of these velocity meas- 
urements is of the order of ±0.5 percent. 

Acceleration chamber shock velocity .- The reported helium shock velocity is the 
average velocity over the distance interval from the secondary diaphragm to the test 
section. The estimated precision of these velocity measurements is of the order of 
±0.2 percent. 

The possibility of a small attenuation of the helium shock velocity was indicated and 
is in agreement with the results reported in reference 1. The shock velocity was also 
measured over a distance interval of about 0.26 m located just upstream of the test sec- 
tion. The estimated precision for these measurements is of the order of ±3 percent. 

The shock velocities obtained in the vicinity of the test section averaged approximately 
2 percent lower than the average velocities over the distance from the secondary dia- 
phragm to the test section. 

Interface velocity .- The reported helium-nitrogen interface velocity was measured 
for a distance interval of about 1.6 m located just upstream of the test section. The esti- 
mated precision of these measurements is of the order of ±0.8 percent. 

An attenuation of the helium -nitrogen interface velocity was not indicated from the 
measurements. This result is in agreement with the results of reference 8, but not with 
the results of reference 1; however, the initial conditions (initial pressures) employed 
here are different than those employed for those two references. Microwave data from 
five tube runs were used to obtain the velocity for 12 equally spaced distance intervals 
located between the secondary diaphragm and the test section. The velocity obtained for 
the interval nearest the diaphragm was about 10 percent lower than the velocities obtained 
for the remaining distance intervals and was expected because of the initial acceleration 
of the interface as reported in reference 1. Except for the distance interval nearest the 
diaphragm, no change in interface velocity was observed to within an uncertainty of 
±1 percent for the remaining distance to the test section. 
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Flow Velocity 


The flow velocity um was obtained from: 

x 2 “ X 1 Ax 
Um = t c> 2 -t c> rAt c 

where x is the position of the center of the photodetector slit in the flow field, and tg 
is the time when the geometric center of the luminous column coincides with the center of 
the photodetector slit. Shown in figure 7 are sample data oscillograms of the two photo- 
detector signal pulses obtained during a velocity measurement. Values of t Cj i and 
t Cj 2 were obtained from these data. 

Model of luminous column .- In order to relate the observed intensity at the photo- 
detectors to spatial position of the luminous column, it was necessary to make at least 
three assumptions concerning the column: the column is symmetrical about its center of 
mass; the total emitted light decays with time; and the column increases in size with 
time. 

As a result of light decay and column growth, the shape of the observed photo- 
detector pulse is asymmetric. The time t p at which peak intensity is observed is 
prior to t c ; that is, 

tp < *c 

Procedure.- The following procedure was used to analyze the data in obtaining 
values of t c . Data from the oscillograms were reduced to numerical form for computa- 
tion at time increments of 0.05 Msec. Corrections were applied to the data for the non- 
linear response of the photodetectors and the decay of the light emitted by the column. 
Sample traces of corrected data pulses are shown in figure 8. Growth of the column can 
be inferred from a comparison of the two pulse widths. 

The remaining asymmetry in the pulse shape was then assumed to be due to the 
growth of the column. The time t c was taken to be the time for which one -half of the 
integrated pulse intensity had been observed. 

Decay correction .- The correction for the decay of the emitted light consisted of the 
reciprocal of an empirically determined decay function. The decay function was derived 
from decay data obtained under static (no flow) conditions for nitrogen densities in the 
range of those encountered in the velocity measurements. The decay data indicated a 
complex decay of the emitted light following excitation of the column. However, the sum 
of two exponential functions approximated the decay during the time interval associated 
with detection of the column at the photodetector slits. The derived decay function used 
in correcting the data was 
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i=0.9e-l-4t + o.le- 0 * 25t 

Io 

where t is the time in psec following column excitation, I is the intensity at time t, 
and Iq is the intensity at time t = 0. 

Approximate size of column .- An initial column diameter of about 1.4 mm based 
on half-peak width was estimated from a photograph of the column obtained under static 
conditions in nitrogen. A short exposure time for the photograph was effected by means 
of a bandpass filter (3920 ± 40 A), which primarily limited the exposure to rapidly decay- 
ing («10 -7 sec) light from the 0-0 band of excited N 2 ions. Shown in figure 9 is a scan of 
the intensity across the luminous column obtained by densitometry from the photograph. 
The location of the scan corresponds to the center line of the expansion tube. Also shown 
for the same location is the estimated diameter of the region illuminated by the ultravio- 
let light source. 

Shown in figure 10 are the estimated diameters of the column as it passed the two 
photodetector slits. These estimates were obtained from the photodetector pulse widths 
and the flow velocity. The time is referenced to the time at which the column was pro- 
duced. A radial velocity of the column boundary of about 150 m/sec can be inferred from 
the results shown in the figure. The observed rate of growth of the column is thought to 
be primarily due to heating of the column gas during excitation. The rate of growth due 
to ambipolar diffusion was estimated to be negligible compared with the observed rate. 

Uncertainties .- Uncertainties associated with Ax and At c lead to uncertainties 
in the flow velocity. The estimated uncertainties in the flow velocity are a random 
uncertainty of about ±1 percent for a single velocity measurement and a systematic 
uncertainty of about ±1 percent for all the velocity measurements. However, it should be 
recognized that the estimate of systematic uncertainty is based mainly on the assumed 
model of the luminous column. The same value of Ax was used in all the measure- 
ments. The estimated uncertainty in - Ax is ±0.5 percent and introduces a systematic 
uncertainty of ±0.5 percent in all the velocity measurements. The nominal value of At c 
was about 4.2 p sec. The estimated uncertainty in the time measurements of approxi- 
mately ±1 percent introduces a random uncertainty in an individual velocity measurement 
of ±1 percent. A systematic error in At c of up to 1 percent may be present as a result 
of errors in the assumed correction procedure. For example, it was assumed that the 
decay correction obtained under static conditions was applicable to the flow conditions. 

The correction procedure was expected to effect only a small correction to estimates of 
Atc which were based on uncorrected data, since the observed photodetector pulse shapes 
were close to being symmetrical. Estimates of At c based on corrected and uncorrected 
data systematically differed by approximately 0.7 percent. 
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Derived Results 


Shocked conditions.- The flow conditions following the passage of the helium and 
nitrogen shocks were derived from the initial conditions and the measured shock veloc- 
ities. The relations used for computation were the adiabatic relations for a normal shock 
given in reference 6. Both helium and nitrogen were assumed to behave as perfect gases. 
The values used for y were 5/3 for helium and 7/5 for nitrogen. 

Density .- The flow density was derived from the measured flow velocity, pitot 
pressure, and tube wall pressure. The values of flow velocity, for times intermediate 
to those for which the velocity was measured, were obtained by linear interpolation 
between the measured values. For times greater than 600 jusec following the arrival of 
the helium shock, velocity values were obtained from a linear extrapolation of the trend 
of the measured velocity. 

The density was computed by using the relation 



M 2 y ’ 

u^ 

_f(M2,yj_ 


where f(M2,y) is the supersonic pitot tube expression given in reference 5. For super- 
sonic flow p « and the term in brackets was considered to be a correction term with 
u^ 

a value close to unity. Approximate values for M were derived from the measured pitot 
pressures and tube wall pressures using the following approximation: 


M 


2 



M lr 

f(M|, r ) 


where 


M 


2 = Vt_ 
a " yiV 


RESULTS AND DISCUSSION 

Average values of the measured and derived results are presented in the figures 
and in tables II and III. The variances presented with the results indicate only the run- 
to-run reproducibility of the results. The symbol a denotes the standard deviation 
from the mean value, and a denotes the estimated variance of the mean value. The 
estimate used was a = u /\} n where n is the number of expansion tube runs from which 
the mean value was obtained. The time is referenced to the start of flow at the test sec- 
tion. In this time frame, the helium shock arrived at t = 0, and the helium -nitrogen 
interface arrived at t = 200 ± 5 jusec. 
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Predicted values from the idealized model are shown in some of the figures for 
comparison with die observed results. Initial conditions used for the model were the 
initial conditions in region @ (fig. 1) and the conditions in region (2) derived from the 
measured nitrogen shock velocity. For computation, equations (1) and (2) and the normal 
shock relations given in reference 6 were used. The arrival time of the nitrogen shock 
at the second diaphragm was taken to be the same as that of the observed flow. For the 
model the helium shock arrived at the test section at t = -360 p,sec and its arrival is 
not indicated in the figures; the helium -nitrogen interface arrived at t = 275 jusec, and 
the expansion fan arrived at t = 72 5 fisec. 

Pitot Pressure 

The observed trend of the pitot pressure is shown in figure 11. During the nitrogen 
flow, departure from a constant pressure is apparent. The large decrease in pressure in 
the interval from about 300 to 450 psec is the pitot pressure dip described in reference 1. 
The minimum pressure of the dip is about 50 percent lower than the preceding peak 
pressure. During the 200 jusec of helium flow, a gradual increase in pitot pressure is 
indicated. The apparent lack of abrupt pressure changes associated with the arrival of 
the helium shock and helium -nitrogen interface is, at least partially, due to the filter- 
ing out of the higher frequency components of the pressure transducer signal. 

The observed trend and magnitude of the pitot pressure are consistent with the 
results presented in reference 1, for which air instead of nitrogen was used as the test 
gas. The trends of pitot pressure, reported in figure 16 of reference 1, for differing 
initial conditions are shown in figure 12. For the initial conditions (figs. 12(a) and 12(b)), 
a pitot pressure dip is apparent, but is not apparent in the trends associated with condi- 
tions of figures 12(c) and 12(d). The initial conditions employed herein most closely cor- 
respond to the conditions shown in figure 12(a). The initial nitrogen pressure employed 
was about twice the initial air pressure employed for conditions of figure 12(a), which 
resulted in pitot pressure levels about twice those shown for conditions of figure 12(a). 
For the same initial pressures of nitrogen or air, the resultant pitot pressure levels are 
expected to be essentially the same. 

Tube Wall Pressure 

The observed trend of the tube wall pressure is shown in figure 13. During the 
nitrogen flow, departure from a constant pressure is apparent, but the relative pressure 
variations are smaller than is the case for the pitot pressure. Some points of similarity 
in the wall pressure and pitot pressure may be noted. However, the wall pressure peak 
and minimum occur earlier, and the minimum pressure is only 20 percent lower than the 
preceding peak pressure. 
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As was noted in reference 1, the wall pressure signal becomes noisy in the vicinity 
of the pitot pressure dip arrival. This can be seen in the sample data oscillogram shown 
in figure 6. The results shown in figure 13 do not reflect the noise oscillations, since 
these high-frequency oscillations were visually averaged in reducing the data oscillograms 
to numerical form. The time interval for which these oscillations occur is also shown in 
figure 13. 

During the 200 psec of helium flow, an increase in wall pressure is indicated. 

This trend is in agreement with the calculations of reference 9. 

The wall pressure predicted by the model is larger than the observed pi'essure for 
the first 580 psec of the real flow. 

Flow Velocity 

The trend of the observed flow velocity is indicated in figure 14. During the nitro- 
gen flow, departure from a constant flow velocity is apparent, but the relative variation 
in velocity is smaller than the relative variations in the observed pressures. All the 
velocities in this time interval can be included within a range of 9 percent. A maximum 
in the observed flow velocity occurs during the passage of the pressure minimum of the 
pitot pressure dip. During the 160-psec interval between the interface and the pressure 
minimum, the velocity increases by 2.5 percent. After the passage of the pressure min- 
imum, the velocity decreases by 9 percent in the next 240 Msec. 

An increase in flow velocity following the arrival of the interface is not unique to 
the operating condition for which a pitot pressure dip is observed. Larger increases in 
velocity (approximately 8 percent) were reported in reference 9 for operating conditions 
which differed significantly from the dip condition. For those conditions the higher veloc- 
ities were observed within a 160-Msec interval following the interface, but locations of 
the maxima were not obtained. The value of pj was 0.5 of the value used here, and the 
values of Pjq were about 0.05, 0.01, and 0.005 of the value employed. 

The helium flow velocity immediately following the helium shock is also shown in 
figure 14. During the 200 Msec of helium flow, the increasing trend of the flow velocity 
(indicated in the figure by short dashed lines) was assumed since no velocity measure- 
ments were made in this interval. The relative values of flow velocity at the contact 
surface (measured) and flow velocity behind the shock (calculated) are in agreement with 
the ratios of the values determined in reference 9. 

The velocity predicted by the model for the nitrogen flow is about 5 percent smaller 
than the velocities observed in the interval from the interface arrival time to 450 Msec. 
The predicted and observed velocities are in agreement late in the flow at 450, 500, and 
600 Msec. The velocity predicted by the model for the helium flow is larger than the 
helium flow velocity immediately following the helium shock. 
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Derived Density 


The density derived from the pitot pressure and flow velocity is shown in figure 15. 
The trend of the density is essentially the same as the trend of the pitot pressure, as 
might be expected from a consideration of the relatively small changes that occur in the 
flow velocity. The dip in pitot pressure corresponds to a dip in density as suggested in 
reference 1 and concluded in reference 8. 

The density predicted from the idealized model is of the order of twice the derived 
density for the nitrogen flow. The predicted arrival of the interface approximately cor- 
responds to the beginning of the pitot pressure dip. 

Also shown in figure 15 is the nitrogen density which is predicted for an isentropic 
expansion from region (2) by using equation (1), the pressure and density in region (2), 
and the observed wall pressure at the test section. The derived density throughout the 
nitrogen flow interval is much less than the density obtained as the result of an isentropic 
process that may occur between region (2 ) and the test section. As inferred from the 
figure, the change in entropy is not constant over the nitrogen flow interval. A trend of 
decreasing entropy change for the flow interval following the density minimum is apparent. 
An extrapolation of this trend suggests that the expansion process may approach an isen- 
tropic process for t > 1000 psec. 

Since the nitrogen expansion process is not isentropic for this expansion tube oper- 
ating condition, any flow models based solely on isentropic processes can not predict the 
derived density at the test section. This result confirms the suggestion made in refer- 
ence 1 that a possible weak (isentropic) P-family expan si on-wave interaction would prob- 
ably not lead to an explanation of the pitot pressure dip or the low values of pitot pressure 
prior to the dip. 

Based on the measurements reported herein, it is not possible to draw any conclu- 
sions concerning the possible mixing of the helium and nitrogen flows as suggested in 
reference 1. In the remaining presentation, it is assumed that no mixing of helium and 
nitrogen occurs for the observed flow. 


Derived Flux 

The timewise integral of the mass flux pu reflects the lower density of the nitro- 
gen flow as compared with the predictions of the idealized model and indicates a relative 
motion between the helium and nitrogen components of the observed flow. For the model, 
the integral from -360 psec to 725 psec is 65.8 g/m^, which consists of 4.2 g/m^ of 
helium and 61.6 g/m2 of nitrogen. For the observed flow the integral from 0 psec to 
715 psec is 37,8 g/m2, which consists of 1.0 g/m^ of helium and 36.8 g/m^ of nitrogen. 
From comparison of the helium component for the two cases, it can be inferred for the 
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observed flow that the nitrogen near the interface has passed through about 76 percent of 
the helium that was initially located between the second diaphragm and the test section. 
The ratio of the helium to nitrogen components for the two cases is 0.068 for the model 
and 0.027 for the observed flow. 


Derived Temperature 

The temperature derived from the observed wall pressure and density by using the 
equation of state is shown in figure 16. A maximum in the observed temperature occurs 
at the same time as the minimum of the pitot pressure dip. The maximum temperature 
of the nitrogen is approximately the same as the temperature of the nitrogen in region ( 2 ). 
The predicted temperature is lower than the observed temperature throughout the nitro- 
gen flow interval. 


Run-to-Run Reproducibility 

The reproducibility of several of the measured flow properties of the expansion 
tube cycle was examined. The measure used in this paper to infer or compare the 
reproducibility is the fractional variance from the mean (that is, o/Mean value). 

Helium- nitrogen interface velocity u^.- Since the interface velocity was used to 
normalize the flow velocity measurements, it was of interest to see whether variations 
in reflected variations in expansion tube operation. It was assumed that variations 
in nitrogen shock velocity U s> i reflected variations in the rupture of the primary dia- 
phragm. For n = 59, the fractional variance of U s> i/uj was ±1.0 percent, with 
corresponding variance for Ug^ i and u^ of ±2.2 percent and 1.9 percent. If the 
variations in Uj were assumed to be unrelated to the variations in U Sj i, the expected 
fractional variance of U s> i/uj would be ±2.9 percent. 

Measured flow velocity u m .- Similar results were obtained for comparisons of the 
variations of uj and the flow velocity u m . For u m corresponding to the peak and 
minimum of the pitot pressure, at t = 275 psec and t = 360 psec, n = 7. For these 
two time locations the fractional variances of u m /ui were ±1.2 percent and ±0.6 per- 
cent with corresponding variances for u m and ui of about ±2.3 percent and 2.2 per- 
cent. If the variations in u m were assumed to be unrelated to the variations in uj, the 
expected fractional variance of u m /ui would be about ±3.2 percent. 

Flow velocity, wall pressure, and pitot pressure u,p w ,p^.- The fractional vari- 
ances '^rthree 7 properHes _ mea^ured _ aFUmTest^^ in figure 17. Listed 

in order from the most to the least reproducible, the properties are (1) the flow veloc- 
ity, (2) the tube wall pressure, and (3) the pitot pressure. 
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Significant increases in the fractional variances occur in the interval of the pitot 
pressure dip for the three properties. The increases occur in the vicinity of the minimum 
of the pitot pressure dip for the flow velocity and pitot pressure, which were observed at 
the center line of the expansion tube, but occurs near the onset of the dip for the tube wall 
pressure. A significant decrease in the fractional variance occurs at approximately 
t = 520 /Li sec for the pitot pressure, but no significant changes are apparent in this 
vicinity for the two other properties. 

Integrated flux .- In practice, the definite integral of a property over a reasonably 
finite time interval is often of more interest than the distribution of that property within 
the time interval, particularly, for an impulse facility such as an expansion tube. The 
fractional variance of the definite integral of the mass flux is shown in figure 18 for sev- 
eral time intervals. The variance of the derived flux pu reflects only the variance of 
the pitot pressure which is shown again for comparison. The relatively smaller vari- 
ance of the velocity was not considered in obtaining the variance of the flux. Although a 
rather large variance for the flux can be expected for an arbitrary time during the nitro- 
gen flow, the integral of the flux for a time interval of 500 psec can be expected to be 
reproducible to within about ±3 percent. 

SUMMARY OF RESULTS 

Measurements of flow velocity, pitot pressure, and tube wall pressure were made 
employing a single expansion tube operating condition, for which a significant decrease 
or dip in pitot pressure occurs during the test flow time interval. Supplemental to the 
three properties measured at the test section were measurements of the intermediate 
chamber shock velocity, the acceleration chamber shock velocity, and the helium- 
nitrogen interface velocity. The following characteristics of the test flow and expansion 
tube flow cycle were derived from these measurements. 

Main Results 

The main results are as follows: 

1. The variations which occur during the test flow time interval in flow velocity 
and tube wall pressure are small relative to the variations in pitot pressure. 

2. In the time interval in which the pitot pressure dip is observed, a corresponding 
dip in the derived flow density occurs. 

3. With regard to the overall expansion processes which occur between the shocked 
conditions of the test gas in the intermediate chamber and the conditions at the test sec- 
tion, a larger expansion occurs than can result from only isentropic processes. 
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Specifically, the value of the average density over the observed nitrogen flow interval is 
about one-half the value that would result from an isentropic expansion to the observed 
pressure. This result is consistent with the lower than expected levels of pitot pressure 
reported in NASA TN D-3421. 


Specific Results 

The specific results are as follows: 

1. The interface velocity measured by means of microwave agrees with the flow 
velocity measured 30 psec following passage of the interface to within an uncertainty of 
about ±1 percent. 

2. Following the passage of the interface, the flow velocity increases to a maximum 
value 2.5 percent larger than the interface velocity during 160 p sec of flow, and then 
decreases to a value about 6.5 percent smaller than the interface velocity during the next 
240 psec. 

3. The maximum in the flow velocity occurs at a time which coincides with the pas- 
sage of the minimum pressure of the pitot pressure dip. 

4. Small variations in the intermediate chamber shock velocity are reflected in the 
interface velocity measured at the test section. 

5. Small variations in the interface velocity are reflected in the nitrogen flow veloc- 
ity in the interval from interface arrival to arrival of the minimum pressure of the pitot 
pressure dip, but may be reflected to a lesser degree for subsequent time intervals. 

6. Of the three properties measured at the test section within the nitrogen flow 
interval, the most reproducible property was the flow velocity; next was the tube wall 
pressure; and the least reproducible was the pitot pressure. 

7. An increase in the fractional variance of the tube wall pressure occurred near the 
onset of the pitot pressure dip, whereas increases in the fractional variances of the pitot 
pressure and the flow velocity occurred nearer the minimum pressure of the pitot 
pressure dip. 

8. The timewise integral of the mass flux over the 515 psec nitrogen flow interval 
was reproducible to within about ±3 percent. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 3, 1974. 
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TABLE I*- INITIAL CONDITIONS 
^Expansion tube diameter, 9.55 cmj 


Chamber 

Length, 

m 

Filling gas 

Pressure, 

kN/m 2 

Temperature, 

K 

Driver 

1.61 

Hydrogen 

8300 

297 

Intermediate 

10.97 

Nitrogen 

6.67 

297 

Acceleration 

9.72 

Helium 

0.267 

297 




















t, p w » 

paec kN/m 2 kN/m 2 

5 7.31 ±0.13 
10 7.54 .12 
15 7.76 .11 
20 7.95 .12 
25 8.16 .13 
30 8.38 .12 
35 8.57 -12 
40 8.77 .12 
45 8.96 .12 
50 9.14 .13 
55 9.28 .13 
60 9.40 .14 
65 9.45 .13 
70 9.43 .12 
75 9.38 .11 
80 9.30 .10 
85 9.21 .09 
90 9.16 .08 
95 9.12 .08 
100 9.09 .09 
105 9.08 .10 
110 9.10 .09 
115 9.14 .09 
120 9.20 .08 
125 9.29 .08 
130 9.42 .09 
135 9.52 .13 
140 9.62 .17 
145 9.64 .16 
150 9.63 .14 
155 9.70 .10 
160 9.75 .08 
165 9.72 .08 
170 9.64 .12 
175 9.50 .13 
180 9.30 .11 
185 9.06 .05 
190 8.89 .05 
195 8.80 .06 
200 8.79 .08 
205 8.83 .09 
210 8.94 .12 
215 9.07 .12 
220 9.19 .11 
225 9.34 .09 
230 9.45 .09 
235 9.52 .10 
240 9.54 .15 
245 9.54 .20 


Tube w 

all press* 

ire; n « 

t. 

V 

0 , 

ysec 

kN/m 2 

kN/m 2 

250 

9.47 

±0.20 

255 

9.18 

.15 

260 

9.08 

.16 

265 

8.96 

.13 

270 

8.77 

.11 

275 

8.52 

.18 

280 

8.26 

.27 

285 

8.06 

.29 

290 

7.83 

.25 

295 

7.75 

.18 

300 

7.83 

.14 

305 

8.15 

.11 

310 

8.43 

.10 

315 

8.65 

.10 

320 

8.74 

.16 

325 

8.65 

.24 

330 

8.58 

.24 

335 

8.50 

.18 

340 

8.49 

.19 

345 

8.37 

.27 

350 

8.27 

.34 

355 

8.25 

.32 

360 

8.23 

.23 

365 

8.26 

.17 

370 

8.38 

.13 ; 

375 

8.51 

.08 

380 

8.62 

.10 

385 

8.69 

.13 

390 

8.79 

.15 

395 

8.89 

.17 

400 

8.97 

.20 

405 

9.05 

.22 

410 

9.12 

I .25 

415 

9.23 

.26 

420 

9.39 

.24 

425 

9.56 

.25 

430 

9.65 

.27 

435 

9.66 

.29 

440 

9.71 

.30 

445 

9.84 

. 26 

450 

9.76 

.26 

455 

9.75 

.30 

460 

9.84 

.22 

465 

9.89 

.10 

470 

9.82 

.09 

475 

9.74 

.14 

480 

9.73 

.21 

485 

9.82 

.22 

490 

9.79 

.27 

495 

9.75 

.32 


P w . 

psec kN/m 2 | kN/m 2 

500 9.80 ±0.26 

505 9.92 .25 

510 10.07 .26 

515 10.20 .30 

520 10.25 .35 

525 10.27 .34 

530 10.31 .24 

535 10.21 .20 

540 10.05 .20 

545 10.03 .20 

550 10.08 .22 

555 10.08 .24 

560 10.17 .27 

565 10.36 .26 

570 10.51 .28 

575 10.70 .30 

580 10.96 .31 

585 11.17 .31 

590 11.20 .33 

595 11.08 .31 

600 10.91 .29 

605 10.77 .27 

610 10.70 .26 

615 10.66 .26 

620 10.64 .28 

625 10.65 .31 

630 10.66 .33 

635 10,73 .27 

640 10.81 .20 

645 10.87 .14 

650 10.94 .36 

655 11.04 .21 

660 11.21 .23 

665 11.39 .25 

670 11.56 .25 

675 11.65 .28 

680 11.69 .31 

685 11.67 .35 

690 11.61 .39 

695 11.59 . .41 

700 11.54 .37 

705 11.45 .31 

710 11.40 .22 

715 11.49 .16 












TABLE II,- MEASURED RESULTS - Concluded 


(c) Boundary velocities; n = 4 


Boundary 

Velocity, 

m/sec 

m/sec 

Nitrogen shock 

2406 

±6 

Helium shock 

4731 

12 

He-N 2 Interface 

4376 

18 


(d) Flow velocity 


t. 

Msec 

u, 

m/sec 

m/sec 

n 

230 

4376 

±20 

3 

275 

4437 

21 

7 

360 

4485 

10 

7 

400 

4293 

60 

3 

450 

4-179 

70 

2 

500 

4197 

77 

2 

600 

4092 

96 

2 





TABLE III,- DERIVED RESULTS 


(a) Conditions following nitrogen shock derived 
from initial conditions and shock velocity 

m/sec 1960 

> m/sec 1110 

P 2 , kN/m 2 364 

P 2 » g/m 3 410 

T 2> K 2990 

(b) Conditions following helium shock derived 
from initial conditions and shock velocity 

u 20 » m/sec 3390 

a 2() , m/sec 2810 

P 2 g, kN/m^ 7.20 

P 20 , g/m 3 1.52 

T 20 , K 2280 
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Figure 1.- Distance-time diagram of an idealized expansion tube flow cycle 

showing pertinent flow regions. 
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Figure 2.- Arrangement of components at expansion tube 
test section. Side view. 
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-HW- Hole diameter 0.78 mm 



Figure 3.- Sketch of windowless ultraviolet light source. 
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Figure 5.- Sketch of optical arrangement. 



Increasing signal + Increasing signal 



(a) Pitot pressure. 



(b) Tube wall pressure. 

Figure 6.- Sample data oscillograms of pitot and tube 
wall pressure gage signals. 



Increasing signals 


o 



TT 

0.5 volt 

V 


IfT 

0 

y 


0.2 volt 


1 ysec 


Figure 7.- Sample data oscillograms of photodetector signals. 




Relative Intensity 





Time, microseconds 

Figure 10.- Estimated diameter of luminous column at the photodetector slit positions. 
Time is referenced to the excitation of the column. 
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re, kN/m 



10 psi 
(69 kN/m 2 ) 


10 psi 


10 psi 


10 psi 



1.6 1.8 2.0 2.2 2.4 2.6 


Time after second diaphragm rupture, milliseconds 


( a ) Pjo = 1.98 torr = 264 N/W*. 

(b) P 10 = 0.95 torr = 127 N/m2. 

(c) PiQ = 0.45 torr = 60.0 N/m2. 

(d) PjQ = 0-095 torr = 12.7 N/m2. 

Figure 12.- Pitot pressure tracings showing effect of varying acceleration-chamber 
pressure p 1Q . (Reproduction of fig. 16, ref. 1.) 
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Figure 14.- Flow velocity. 


40 






wv 


Helium flow 

Nitrogen flow 

Model 


I i l I I ! L 

0 200 400 600 

Time, microseconds 


Figure 16.- Temperature. 
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Figure 18.- Variance of pitot pressure and definite integral of mass flux. 
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